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■ Abstract 

The exclusive photoproduction of J/ip mesons, ■yp J/'ipp, has been stud- 
^ I ied in ep coUisions with the ZEUS detector at HERA, in the kinematic range 

O I 20 < W < 290 GeV, where W is the photon-proton centre-of-mass energy. The 

>< I J/ip mesons were reconstructed in the muon and the electron decay channels 

V I using integrated luminosities of 38pb~^ and 55pb~^, respectively. The helicity 

structure of J/ production shows that the hypothesis of s-channel helicity con- 
' servation is satisfied within two standard deviations. The total cross section and 

the differential cross-section da/dt, where t is the squared four-momentum trans- 
^ I fer at the proton vertex, are presented as a function of W, for \t\ < 1.8 GeV^. The 

5^ \ t distribution exhibits an exponential shape with a slope parameter increasing 

logarithmically with W with a value b = 4.15 ± 0. 05 {stat.)tQil{syst.) GeV"^ at 
W = 90 GeV. The effective parameters of the Pomeron trajectory are aip(O) = 
1.200±0.009{stat.)toToisyst.) and a'jp = 0.115±O.Ol8{stat.)to^[l{syst.) GeV~l 
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1 Introduction 



Exclusive J/ip photoproduction is expected to be described by models based on pertur- 
bative QCD (pQCD), since the mass of the charm quark provides a hard scale [1-9]. In 
such models, the photon fluctuates into a cc pair which subsequently interacts with the 
proton. This interaction is modelled by the exchange of a gluon ladder and the cross 
section is proportional to the square of the gluon density. These models predict a rapid 
rise in the cross section with IV, where W is the photon-proton centre-of-mass energy, 
which is caused by the fast increase of the gluon density in the proton at the small values 
of Bjorken x. 

Within the framework of Regge phenomenology [10], diffractive interactions at large 
centre-of-mass energies are the result of the t-channel exchange of the Pomeron trajectory, 
Q!jp(t), carrying the quantum numbers of the vacuum. The differential cross section at 
high energies is expressed as 

^ocF(t)-iy^[«^W-i], (1) 

where t is the squared four-momentum transfer at the proton vertex and F{t) is a function 
of t only. If da/dt decreases exponentially and the trajectory is linear in t, apit) = 
ttip(O) -|- a'jpt, the cross section can be expressed as 

— cx . e''^^^*, (2) 

dt 

where the slope parameter is b{W) = &o + 4:a'jp\n(W /Wq) and Wq is an arbitrary energy 
scale parameter. A flt to hadronic data [11] yields the soft-Pomeron parameter a'jp = 
0.25 GeV~^. The pQCD models predict the effective a'jp in the perturbative regime [2,6, 
12, 13] to be much smaller than 0.25 GeV^. 

Studies of the exclusive, diffractive photoproduction of vector charmonium states, i.e. 
7P J/'ipP, at HERA show that the total cross section [14,15] rises steeply with W. 
In addition, the t dependence of the cross section can be fltted by a single exponential, 
da/dt ~ e~^'*', with b ~ 4.6 GeV~^ [14,15]. There are indications that the slope parameter 
b has little variation with W [15,16], i.e. little "shrinkage" is observed. 

In this paper, more precise measurements of the t and W dependence of the exclusive 
photoproduction cross section J/ip mesons are made using the reaction ep — e7p — >■ 
eJ/ipp for values of the photon virtuality, Q^, close to zero in the range 20 < W < 
290 GeV. With respect to the previous ZEUS result [14], this analysis covers wider ranges 
in W and t and has a large increase in statistics, combined with an improved understanding 
of the detector and of the background subtraction. The cross section for 7p — s> J/ipp 
and the slope of the differential cross-section da^p^j/^p/dt are studied as a function of W. 
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The Pomeron trajectory parameters ajp{0) and a'jp, which describe the W dependence of 
the cross section and the shrinkage, respectively, are determined. The hehcity structure 
of J/ip production is investigated to test the vahdity of SCHC. 

2 Experimental set-up 

In this analysis, J/ip mesons were identified with the ZEUS detector at HERA by their 
decays to or to e~^e~. The muon sample corresponds to an integrated luminosity 

of 38.0 ± 0.6 pb~^, collected in 1996 and 1997 when 27.5 GeV positrons^ were collided 
with 820 GeV protons. For the measurement with 20 < < 30 GeV, a sample from an 
integrated luminosity of 27.5 ± 0.4 pb~^ was used. The electron sample corresponds to 
an integrated luminosity of 55.2 ± 1.2 pb~^, collected in 1999 and 2000 in collisions of 
27.5 GeV positrons with 920 GeV protons. 

The ZEUS detector is described in detail elsewhere [17]; only the components most rele- 
vant for this analysis are outlined here. 

Charged particles are tracked by the central tracking detector (CTD) [18], which operates 
in a magnetic field of 1.43 T provided by a thin super-conducting coil. The CTD consists 
of 72 cylindrical drift chamber layers, organised in 9 superlayers covering the polar-angle^ 
region 15° < 9 < 164°. The relative transverse-momentum resolution for full-length tracks 
is a{pT)/pT = 0.0058pT © 0.0065 © 0.0014/pt, with pt in GeV. 

Charged particles in the forward direction are detected in the forward tracking detector 
(FTD) [17], which consists of three planar drift chambers perpendicular to the beam 
covering the polar angles 7° < 6 < 28°. Each chamber is made of three layers of drift 
cells; the wire directions in the second layer are rotated by 120° with respect to the first 
layer and similarly for the third layer with respect to the second. Each drift cell has six 
sense wires. Thus each chamber measures a track segment in three projections with up 
to six hits per projection. 

Surrounding the solenoid is the high- resolution uranium-scintillator calorimeter (GAL) [19]. 
It consists of three parts: the forward (FCAL, 2.6° < 9 < 36.7°), the barrel (BCAL, 

^ Hereafter, "positron" is used to refer to both electron and positron beams. At the values of 
studied here, and e'^p scattering were assumed to give identical results since contributions from 
Z'^ exchange are negligible. Similarly, "electron" is used to refer to either the electron or positron from 
the decay of the J/ip- 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z-axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X-axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as 77 = — ln(tan |), where the polar angle, d, is measured with respect to the proton beam 
direction. 
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36.7° < 9 < 129.1°) and the rear (RCAL, 129.1° < 9 < 176.2°) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). 
The smallest subdivision of the calorimeter is called a cell. Under test beam conditions, 
the CAL has relative energy resolutions of a{E)/E = 0.18/^/E for electrons hitting the 
center of a calorimeter cell and a{E)/E = 0.35 /VE for single hadrons {E in GeV). Cell 
clusters were used to aid in the identification of muons and electrons. 

The forward plug calorimeter (FPC) [20] is a lead-scintillator sandwich calorimeter with 
wavelength-shifter fibre readout. Installed in 1998 in the 20 x 20 cm^ beamhole of the 
FCAL, it has a small hole of radius 3.15 cm in the centre to accommodate the beampipe. 
It extends the pseudorapidity coverage of the forward calorimeter from 77 < 4.0 to 77 < 5.0. 

The small-angle rear tracking detector (SRTD) [21] is attached to the front face of the 
RCAL. The SRTD consists of two planes of scintillator strips read out via optical fibres 
and photomultiplier tubes. It covers the region 68 x 68 cm^ in X and Y with the exclusion 
of a 8 X 20 cm^ hole at the centre for the beampipe. The SRTD provides a transverse 
position resolution of 3 mm and was used to measure the positions of electrons, from the 
J/ ip decay, produced at small angles to the positron beam direction. 

The muon system consists of limited streamer tubes (forward, barrel and rear muon 
chambers: FMUON [17], B/RMUON [22]) placed inside and outside the magnet yoke. 
The inner chambers, F/B/RMUI, cover the polar angles between 10° < 9 < 35°, 34° < 
9 < 135° and 135° <9 < 171°, respectively. The FMUON has additional drift chambers 
and permits high-momentum muon reconstruction for polar angles between 6° and 30° 
using the magnetic field of 1.7 T produced by two iron toroids placed at Z = 9 m and 
the toroidal field of 1.6 T provided by the yoke coils. The relative momentum resolution 
of o'{p)/p = 0.2, up to 20 GeV, is dominated by the multiple scattering. 

The proton-remnant tagger (PRTl) [23] consists of two layers of scintillation counters 
located at Z = 5.15 m, and covers the pseudorapidity range 4.3 < 77 < 5.8. It was used, 
up to the end of the 1997 running period, to tag events in which the proton diffractively 
dissociates. 

The luminosity was determined from the rate of the bremsstrahlung process ep — > e'jp, 
where the photon was measured with a lead-scintillator calorimeter [24] located at Z = 
-107 m. 
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3 Kinematic variables and reconstruction 



The kinematic variables used to describe exclusive J/ip production, 

e{k)p{P) e{k') J/tP{v)p{P'), 

where k, k', P, P' and v are the four-momenta of the incident positron, scattered positron, 
incident proton, scattered proton and J/ip meson, respectively, are defined as follows: 

• Q"^ = —q^ = —[k—k'y, the negative four-momentum squared of the exchanged photon; 

• = {q + Py, the squared centre-of-mass energy of the photon-proton system; 

• y = {P ■ q)/{P ■ k), the fraction of the positron energy transferred to the photon in the 
rest frame of the initial-state proton; 

• t = {P — P'Y = — the square of the four- momentum transfer at the proton 
vertex. 

The following angles are used to describe the decay of the J/ip: 

• Oh and (phi the polar and azimuthal angles of the positively charged decay lepton in 
the helicity frame, defined as the J /if) rest frame with the quantisation axis taken to 
be the J/ip direction in the photon-proton centre-of-mass system. The origin of the 
azimuthal angle, (ph, corresponds to the case when the decay particles are produced in 
the production plane, defined as the plane containing the incident photon (assumed 
to be in the positron direction) and the J/ip momentum vectors. 

In addition to the above quantities. My, the mass of the diffractively produced state 
y, characterises the major background process, the proton-dissociative reaction ep — > 
e J/i) Y. 

For the selected events, ranged from the kinematic minimum, = Mly^ j{ \ —y)^, 

10"^^ GeV^, where Mg is the positron mass, up to Ql^ax ~ ^ GeV^, the value at which the 
scattered positron starts to be observed in the GAL, with a median of approximately 
5 ■ 10"^ GeV^. Since the typical is small, it can be neglected in the reconstruction of 
the other kinematic variables. 

The photon-proton centre-of-mass energy, W , can be expressed as ~ 2Ep(E — pz)ii, 
where Ep is the laboratory energy of the incoming proton and {E —pz)ii is the difference 
between the energy and the longitudinal momentum of the dilepton system. 

The squared four-momentum transfer at the proton vertex is given by t ^ ~Pt^ the 
transverse-momentum squared of the dilepton system. Non-zero values of give t values 
that differ from — by less than Q^; this effect is corrected for using the Monte Garlo 
simulation. 
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Since neither the scattered positron nor the scattered proton was observed, the kinematic 
variables were reconstructed using only the measured momenta of the decay particles. At 
low W, the J/ip mesons are produced in the forward direction, while at high W they are 
produced in the backward direction. For the J/ip ^ fi~^fi~ sample, the CTD, FTD and 
FMUON information were used when available. For the J/tp ^ e~^e~ sample, the CTD 
or the CAL/SRTD information was used. The relevant in situ electron-energy resolution 
of the CAL for energies in the range 3 < -E < 20 GeV averages a{E)/E = 0.27 /^/E for 
this analysis [25]. 

4 Trigger and event selection 

The events were selected online via a three-level trigger system [17]. The signature for 
exclusive J/ ip photoproduction events consists of a pair of charged leptons, with no other 
significant activity in either the CTD or the CAL, since the scattered positron and proton 
escape undetected down the beampipe at small scattering angles. For the two decay 
channels, the elasticity cuts described below were imposed. They restrict the photon 
virtuality to Q'^ < 1 GeV^ and the mass of the dissociative system to My ^ 3.0 GeV. 
No cut on t was applied on either channel. To select candidate events for the muon and 
electron decay channels, different selection cuts were applied. 

4.1 Muon channel 

The trigger selected events with at least two tracks in the CTD or one track in the CTD 
and one track in FMUON. At least one track had to point towards an energy deposit 
compatible with a minimum ionising particle (m.i.p.) in the CAL and either a hit in the 
FMUI or a segment in the B/RMUI. The trigger efficiency for events that passed the 
offline selection cuts, defined below, was determined from independent triggers and from 
MC simulations to be between 50% and 75%, depending on W. 

Events having the following characteristics were selected offline: 

• exactly two oppositely charged tracks from a common vertex, with Z-coordinate 
\Zvertex\ < 50 cm, at least one of which matches either a hit in the FMUI or a segment 
in the B/RMUI; 

• each CTD track passes through at least three superlayers, effectively limiting the 
polar- angle region of these "CTD tracks" to 17° < ^ < 163°; 

• the angle between the two tracks is less than 174.2°, in order to reject cosmic-ray 
events; 
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• CAL energy associated with each track consistent with the energy deposit of a m.i.p., 
i.e. between 0.8 and 5 GeV, with a ratio of at least 0.8 between the energies in the 
HAG and the EMC sections. The energy was associated with the track if it was inside 
a cone of radius 30 cm (in EMC) or 50 cm (in HAC), centred at the impact position 
of the track extrapolated on to the CAL; 

• no CAL cell, apart from those associated with a candidate muon, with energy above 
the threshold level of 150 MeV to 200 MeV, depending on the calorimeter part and 
section. This elasticity requirement rejects proton-dissociative and inelastic events as 
well as deep inelastic scattering (DIS) events. The PRTl was not used to veto the 
events. 

To define the kinematic region in which both decay muons could be well reconstructed, the 
analysis was limited to the range 20 < < 170 GeV. The requirement that CTD tracks 
traverse at least three superlayers leads to a small acceptance for two such tracks in the 
region 20 < < 30 GeV. In this region, therefore, additional events were accepted with 
one CTD track and a second track in the FMUGN spectrometer, which both measures 
the momentum and triggers the event. The FMUON track was re-fitted to the vertex 
including the FTD segments to improve the reconstruction parameters. In the kinematic 
region 30 < TV < 170 GeV, only events with two CTD tracks were used. 

4.2 Electron channel 

The J/if) ^ e^e~ events were selected online using two different trigger algorithms: 

• the first algorithm was optimised for events with one or two tracks in the CTD. It 
required at least one track, but fewer than five tracks, in the CTD, a total energy 
deposit greater than 1.8 GeV in an EMC section of the CAL and an energy of less 
than 3.75 GeV in the region of 60 x 60 cm^ of FCAL around the forward beam pipe; 

• the second algorithm was optimised for events with zero tracks in the CTD. It required 
at least two deposits in the RCAL EMC, each with an energy greater than 2.1 GeV. 

The trigger efficiency for events that passed the offline selection cuts, defined below, was 
determined from independent triggers and from MC simulations to be between 80% and 
90%, depending on W . 

The following offline selection requirements were applied: 

• for events with a reconstructed tracking vertex, the cut \ZyTx\ < 50 cm was ap- 
plied. Events without a vertex were accepted and were assigned a vertex Z position 
corresponding to the nominal interaction point; 
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• the CAL energy deposits were grouped into clusters. Events were selected for further 
analysis if they satisfied one of the following criteria: 

— for events with two CTD tracks, the CAL energy associated with each track, inside 
a cone of radius 25 cm centred at the impact position of the track extrapolated 
on to the CAL, had to be consistent with the energy deposition expected for 
an electron: Eemc/E > 0.9, where E and Eemc are, respectively, the total 
energy and the energy deposited in the EMC section. These two tracks were then 
considered electron candidates; 

— for events with one CTD track and one cluster not associated with it, the CAL 
energy associated with the track had to have Eemc/E > 0.9 and the CAL energy 
of the unassociated cluster had to be more than 3 GeV, with Eemc/E > 0.98, to 
be considered an electron candidate; 

— for events with no CTD track, two clusters, each with CAL energy more than 
3.5 GeV and with Eemc/E > 0.98, were required as candidate electrons; 

• any energy deposit in the CAL cells, not associated with either of the two electron 
candidates, was required to be less than either 200 MeV or 300 MeV, depending on 
the calorimeter part and section. This requirement^ rejects proton-dissociative and 
inelastic events, as well as DIS events; 

• further to reduce events from proton dissociation, the energy measured in the FPC 
was required to be less than 1 GeV. 

The analysis was restricted to the kinematic region 20 < < 290 GeV because the 
acceptance drops at lower W and the QED-Compton background dominates at higher 
W. 

5 Monte Carlo simulation 

The acceptance and the effects of the detector response were determined using samples 
of Monte Carlo (MC) events. The ZEUS detector response to the generated particles was 
simulated in detail using a program based on GEANT3.13 [26]. All the generated events 
were processed through the same reconstruction and analysis chain as the data. 

The exclusive processes ep — eJ/ipp and ep eip{2S)p were modelled using the 
MC generators DIPSI [1,27] and ZEUSVM [28]. The events were weighted with a jp 
cross section proportional to and with an exponential t dependence e''*. The weight 

^ This elasticity requirement is less stringent than the one used for the selection of the muon sample 
because of an increased noise level in the CAL during the data taking in 1999 and 2000. 
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parameters, 5 = 0.70 and b = 4.3 GeV~^, were chosen so as to describe the W and 
t dependence of the data, as discussed in Sections 9 and 10. The decay of the J/ip 
mesons in the centre-of-mass system was generated with a (1 + cos^ 6^) distribution, 
consistent with the measurement presented in Section 8. The effect of the initial-state 
radiation on the acceptance, estimated using HERACLES 4.6.1 [29,30], was 5% in the 
region 20 < TV < 30 GeV and neghgible at larger W. The final-state radiation of hard 
photons at large angles with respect to either of the two decay leptons decreases the 
acceptance, as these events are rejected by the elasticity requirements applied offline. 
The effect was estimated using PHOTOS [30,31] and was found to be 1.3% in the muon 
analysis and up to 8% in the electron analysis. The acceptances. A, corrected for these 
effects, are given in Table 1. They were calculated as the number of events reconstructed 
in a bin divided by the number of events generated in the same bin. 

Proton-dissociative events, ep — > eJ/ipY, were modelled using the generator EPSOFT 
[32], which simulates 7p interactions assuming the exchange of the soft-Pomeron trajec- 
tory. At fixed W and t, it models the mass spectrum of the baryon system Y according 
to da/dMy oc My^ , with My > 1.25 GeV. The multiplicity distribution of hadrons from 
the decay of the proton-dissociative system and their transverse momenta with respect 
to the proton-Pomeron collisions axis are simulated to describe the ZEUS photoproduc- 
tion data and hadron-hadron single-diffractive results. The longitudinal momenta are 
generated using a uniform distribution in rapidity. The simulation parameters 5 = 0.70, 
b = 0.65GeV-2 and (3 = 2.6 were chosen to describe the W, t and My dependence 
of the data, as discussed in Section 6.3. For the region My < 2 GeV, b = 4.0 GeV~^ 
was used. This reflects the steeper t distribution observed in low-mass hadron-hadron 
diffraction [33,34]. 

The QED 77 background was simulated using the LPAIR [35] generator. The 

QED-Compton scattering, ep — e'-fp, background was simulated using the COMP- 
T0N2 [36] generator. 

6 Mass spectra and background subtraction 

The invariant-mass spectra for the muon- and electron-pair candidates, after all offline 
cuts, are shown in Figs. 1 and 2 for representative W bins. The mass resolution is excellent 
in the region 30 < < 150 GeV, where both reconstructed leptons pass through all CTD 
layers, and decreases at lower and higher W values, where the leptons are produced in 
the forward and rear direction, respectively, at the edges of the CTD acceptance. In the 
analysis of the electron decay channel, the kinematic region was extended at very high W 
by reconstructing the electrons in the RCAL and SRTD. The mass spectra of the electron 
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pairs have a tail at low mass due to photon bremsstrahlung. 

The final samples contain backgrounds from non-resonant sources: QED 77 processes, 
QED-Compton scattering and other non-resonant background (dominated by misidenti- 
fied pion production), as well as from resonant processes: diffractive iIj{2S) production and 
proton dissociation. Contamination from non-diffractive J/^ production, ep — eJ/ipX, 
estimated using the MC generator HERWIG 5.8 [37], is negligible. 

6.1 Non-resonant background 

Background from the QED process 77 l^l', in which a lepton pair is produced by 
the fusion of a photon radiated by the positron with a photon radiated by the proton, 
was estimated using the LPAIR MC normalised to the data in the mass region outside 
the J/ip and ip{2S) resonances. It is shown in the dilepton-mass spectra of Figs. 1 and 2. 
The typical contribution in the signal region is 10% for the muon sample and up to 20% 
for the electron sample. 

Pions produced at low angles in the forward direction can be misidentified as muons or 
electrons. This background dominates at low masses and for W < 50 GeV, as can be seen 
in Figs. 1(a) and 2(a,b). 

For the electron sample, there is an additional contribution from QED-Compton scattering 
with initial-state radiation. It was estimated using the C0MPT0N2 MC normalised to 
the data in the mass region outside the resonances for W > 230 GeV, where it dominates 
the background distribution, as shown in Fig. 2. The size of this background ranges from 
3% at W = 200 GeV to 50% at W = 275 GeV. The total non-resonant contributions are 
given in Table 1. 

6.2 Events from ^(25) production 

Events from ip(2S) production can fake exclusive J/ip events, mainly through the decay 
il^{2S) J/ip + neutrals, with a branching ratio E = (23.1 ± 2.3)% [38], when the 
J/ip decays into two leptons and the neutral particles are not detected in the CAL. 
An additional source of background comes directly from tp{2S) decays, with a 

branching ratio B = (1.03 ± 0.35)% [38], since the mass window used to count the J/ip 
signal, for W ^ 35 GeV and W ^ 140 GeV, is large enough to include the ^^{23) mass 
at 3.685 GeV. The number of events from 'ip{2S) production present in the J/ip elastic 
sample was estimated using the ip{2S) DIPSI MC events and the ratio of production cross 
sections, ip{2S)/{J/^) = 0.150 ± 0.027{stat.) ± 0.022{syst.) [39], to be smaller than 7%, 
as shown in Table 1. 
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6.3 Proton dissociation 



The largest source of background is given by the diffractive production of J/ if) mesons 
with proton dissociation, ep — > eJ/ipY, when the system Y has a small mass and its 
decay products are not detected in either the FCAL, the PRTl or the FPC 

To estimate this background, the elasticity cut was removed in the region of fragmentation 
of the system Y\ the proton-dissociative data obtained in this way were used to tune the 
EPSOFT MC generator. The data and the fraction of MC events in which energy was 
deposited in the PRTl or the FPC were then used to estimate the number of proton- 
dissociative events contaminating the exclusive J / ip sample. 

The parameters h and (3 of the EPSOFT MC generator were tuned using a sample of 
dimuon events triggered by the B/RMUI. These data were selected as described in Section 
4.1 but, for those FCAL cells with 9 < 30°, the elasticity requirement was removed 
and an energy of at least 300 MeV in both EMC and HAC sections was required. A 
rapidity gap of at least Arj > 1.3 between the J/ip meson and the products of the proton- 
dissociative system was required. The final sample of J ftp candidates, in the kinematic 
region 90 < < 130 GeV, consisted of about 600 events. The data sample corresponds, 
according to MC simulations, to the region 3.5 < My < 30 GeV and p'^ < 10 GeV^. The b 
slope and P parameters were determined to be 0.65±0.10 GeV~^ and 2.6±0.3, respectively, 
from the study of the p^ and visible CAL-energy distributions. The W distribution for the 
proton-dissociative process was consistent with that of the exclusive channel. A similar 
study performed for the electron-decay channel, using a smaller sample of data triggered 
without the veto on energy deposited in the FCAL region around the beam pipe, yielded 
values of b and (3 consistent with those determined from the muon sample. 

The proton-dissociative events misidentified as exclusive J/ ip production were subtracted 
in each {W, t) bin for the cross sections presented here, for both the muon and the electron 
analyses. 

In the muon analysis, the contribution from proton-dissociative events in the elastic 
sample was estimated using the PRTl. In each W and t bin, the quantity fp-diss = 
fpRTi ■ "^^s computed, where /pr^i is the fraction of the data tagged in the PRTl 
and e is the MC tagging efficiency, defined as the probability to obtain a tag in PRTl 
in EPSOFT events. Using /^'^^^ = 12.6% and e = 57.2%, fp.diss was estimated to be 
{22.0 ± 2. {stat.)± 2. 0{sy St. ))% for |t| < 2 GeV^ with no W dependence. The fraction of 
fp-diss was estimated to increase with t from {11.0'ti l{stat. + syst.))%, in the first t bin, 
up to (49+^^(stat. + syst.))%, for 1.2 < |t| < 1.8 GeV^ 

In 1998, the PRTl was no longer used. The FPC was inserted and was used to veto proton- 
dissociative events. In the electron analysis, which uses the data taken in 1999 and 2000, 
the amount of proton-dissociative background in the elastic sample was estimated from 
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fp-diss = f^tc ■ (7 - 1) and found to be (17.5 ±l2{stat.)tll{syst.))%, using fp^^^ = 12% 
for the FPC-tagged events and e = 40.7%, where e is defined as the probabihty to obtain 
a tag in FPC in EPSOFT events. The formula is different with respect to the muon case 
because the events tagged by the FPC were also rejected. The fraction of fp-diss was 
estimated to increase with t from {'6'^\{stat. + syst.))% at low t to {2S't^{stat. + syst.))% 
at 0.85 < t < 1.15 GeVl 

As an independent check on the estimation of the proton-dissociative background, an 
alternative model to EPSOFT was used. In the baryon resonance region, at low My, 
a resonant component with slope h = 6.5 GeV~^ was considered. A second component 
due to non-resonant proton dissociation with slope b = 0.65 GeV~^ was added. The two 
components were constrained to satisfy the first moment of the finite- mass sum rule [40]. 
This model yielded results that agreed with those from EPSOFT to within 2% for both 
the cross section and the b slope. 



6.4 Signal determination 

Since the mass spectra have shapes and background contributions that vary with W, t 
and the decay channel, different procedures were used in the muon and electron analyses 
to determine the number of signal events. 

For the muon analysis, the signal events were counted in a mass window corresponding 
to ± three standard deviations of the Gaussian fit from the mean fitted value of the J/ ip 
mass. Since the mass resolution depends on the kinematic region, the mass windows 
were different for different W and t bins. The number of 77 — fi~^fi~ background 
events was estimated as described in Section 6.1 and subtracted in each mass window to 
obtain the number of J/ip candidate events. For the lowest W bin, there is a remaining 
background, coming from misidentified pions. The mass spectrum was then fitted to the 
sum of the distribution predicted by the signal MC and a single exponential function for 
the background. 

For the electron analysis, the backgrounds due to QED-Compton scattering and from 
77 e^e~ were subtracted. The remaining mass spectrum was then fitted, in each 
bin, to the sum of the distribution predicted by the signal MC and a single exponential 
function for the remaining background. The latter is dominated by misidentified pions. 
The MC gives a good description of the Mg+g- distributions observed in the data with 
typical x^/ndf better than 1.5 
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7 Systematic uncertainties 



The systematic uncertainties on the cross sections are given separately for the two decay 
channels. For the muon channel, the following sources of uncertainty were considered. 

• trigger efficiency: the uncertainty due to that of the trigger efficiency was ±5% for the 
CTD track reconstruction and up to ±7% for the muon selection; 

• event selection: the minimum number of required CTD superlayers was raised from 
three to four; the cut on the angle between two tracks was relaxed from 174.2° to 176.4°; 
the criteria for the identification of a m.i.p. in the CAL were varied and the cell energy 
threshold for the selection of exclusive events was increased from 150 or 200 MeV to 300 
MeV. The resulting uncertainty was ±3%. For events in the range 20 < < 30 GeV, 
the uncertainty due to the track reconstruction in the FMUON/FTD was ±11%; 

• MC model dependence: the uncertainty was estimated by varying the parameters 
b and S of the DIPSI MC simulation within the range 4.1 < b < 4.5GeV~^ and 
0.60 < 6 < 0.75. The centre-of-mass decay angular distribution of the muons was 
changed to [1 + a + (1 — 3a) cos^ 6h] with a = —0.05, consistent with the measurement 
presented in Section 8. The overall uncertainty due to model dependence was ±5%; 

• proton-dissociative subtraction: the uncertainty on the modelling of the hadronic final 
state in proton-dissociative events was estimated to be ±2% by varying the parameters 
of the simulation by ±0.3 for p and ±0.10 GeV~^ (±2 GeV"^) for b when My > 2 GeV 
(My < 2 GeV), as discussed in Section 6; 

• non-resonant background subtraction: the uncertainty was typically ±2%. In the 
lowest W bin, the uncertainty was 7%, as determined by the fit; ; 

• the uncertainty in the luminosity determination was ±1.7% for the 1996-1997 running 
period. 

For the electron channel, the following sources of uncertainty were considered: 

• trigger efficiency: the estimated uncertainty was ±2.5% for the CTD, ±(1 — 5)%, de- 
pending on W, for the CAL energy threshold and ±3% for the trigger stream requiring 
two isolated electromagnetic clusters; 

• event selection: the effect of varying the elasticity requirements by raising the threshold 
by 100 MeV was -2.5% to ±2.5%, depending on W; 

• MC model dependence: the uncertainty, estimated in the same way as described for 
the muon channel, was less than ±2.5%; 

• proton-dissociative subtraction: the uncertainty, estimated as described for the muon 
channel, was 
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• non-resonant background subtraction: the uncertainty in the normahsation of this 
background, as determined by the fit, varied between 1 and 6%, depending on W; 

• the uncertainty in the luminosity determination was ±2.25% in the years 1999 and 
2000. 

The overall systematic uncertainty was determined by adding the uncertainties in quadra- 
ture. An additional uncertainty of 1.7% [38] associated with the branching ratio Bj/^^i+i- 
was not included. Since the major sources of systematic uncertainty are mostly indepen- 
dent of W , they have a small influence on the determination of 5 and the Pomeron 
trajectory. 



8 Decay angular distributions 

Since the decay angular distributions were used to reweight the MC simulated events and 
thus affect the cross-section measurements, they are discussed first. They were used to 
investigate the helicity structure of J /if) production. The decay angular distribution is 
a function of 9h and the polar and azimuthal angles of the positively charged lepton 
in the helicity frame. The normalised angular distributions can be expressed [41] in the 
form 

' \[l+rZ+{l-Kt)co.^9,] (3) 



N d cos Bh 8 
and 

1 dN 1 



d(hh 2ti 



[l + r°^cos20;,] , (4) 



where the J/ip spin-density matrix element Tqq represents the probability that the pro- 
duced J /if) has helicity and r^l^^ is related to the interference between the non-flip and 
double-flip amplitudes. If the J/ip retains the helicity of the almost-real photon, as in 
the hypothesis of s-channel helicity conservation (SCHC), Tqq and r\^^ should both be 
approximately zero. 

The angular distributions of the leptons from J/ip decay are presented in Fig. 3. They 
were measured in the kinematic range 30 < < 170 GeV and |t| < 1 GeV^, using events 
in which both leptons were measured in the CTD. The non-resonant background was 
subtracted in each angular bin. No subtraction of the dissociative contribution was made, 
since the proton-dissociative sample, discussed in Section 6.3, displayed similar angular 
distributions to the elastic events. The ip(2S) events were assumed to have the same 
angular distribution as the J/ip events and were not subtracted. 



13 



The elements Tqq and r^^i, obtained by fitting the acceptance-corrected 9h and (j)h distri- 
butions to Eqs. (3) and (4), are 

= -0.017 ± 0.015(stat.) ± 0.009{syst.) 

and 

r^i^ = -0.027 ± O.Q13{stat.) ± 0.005(syst.). 

Thus, to within two standard deviations, the SCHC hypothesis holds, as expected for 
heavy mesons [42]. 



9 W dependence of the cross section 

The 7p cross section for exclusive J/ ip production was evaluated from the ep cross section 
using the expression 

_ (^ep^ej/i^p _ 1 {Nobs — Nnon-res — N^(-2S)) ' (1 " fp~diss) 

^.v^Jl^v - ^ - ^ JTa^ ' 

where $ is the effective photon flux [43] , No^s is the number of events in the signal mass 
region, Nnon-res is the number of non-resonant background events, N^(2S) is the number 
of events from ip{2S) production, fp-diss is the fraction of proton-dissociative events, £ 
is the integrated luminosity, A is the acceptance and B is the branching ratio, where 
B = (5.93±0.10)% for the electron channel and (5.88±0.10)% for the muon channel [38]. 

The numbers of events, the acceptance, the flux factors and the cross sections are given in 
W bins for each decay mode in Table 1 . The cross section is shown as a function of W in 
Fig. 4. No cut on t was applied. The small difference (~ 4%) in the normalisation of the 
muon and the electron values is within the correlated uncertainty associated with each 
decay channel. The values are larger than those determined previously by ZEUS [14]. 
The differences are due to a better understanding of the acceptance and trigger efficiency 
and of the background subtraction. Because of these improvements, the results of this 
paper supersede those of the previous publication. 

Results from the HI Collaboration [15] and from fixed-target experiments [44,45] are also 
displayed in Fig. 4. While the W dependence is similar, there is a normalisation difference 
between the HI and ZEUS values. 

The results of fits of the form cr oc (1^/90 GeV)^ to the muon data and, separately, to 
the electron data are given in Table 2. A common fit to the data with W > 30 GeV, 
including both the muon and the electron measurements, with statistical and systematic 
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uncertainties added in quadrature, yields a value of 6 = 0.69 ± 0.02{stat.) ± 0.03{syst.). 
This result, shown as the curve in Fig. 4, confirms the strong energy dependence of the 
cross section observed previously [14,15]. The measurements for W < 30 GeV were not 
included in this fit to avoid possible effects due to the charm production threshold [46]. 
However, the result of a fit including the points with W < 30 GeV does not significantly 
change the fitted value of 6. 

The ZEUS data are compared in Fig. 5 to leading-log-approximation (LLA) pQCD calcula- 
tions [5] , based on open cc production and parton-hadron duality, using the CTEQ5M [47] 
(dashed curve) or MRST99 [48] (dotted curve) parton-density functions (PDF). The gluon 
density is evolved using "skewed" evolution equations [49]. 

The solid curve in Fig. 5 is the result of a LLA pQCD calculation [6] based on the 
interaction of the proton with qq dipoles with small transverse size via two-gluon exchange. 
The model uses the CTEQ4L gluon PDF [50] evolved using skewed evolution. This 
calculation is sensitive to the value of A, a scaling parameter that relates the transverse 
size of the dipole to the four-momentum scales in the interaction cross section. The curve 
shown uses A =4, which gives a dependence on W that is less steep than for A = 10, which 
was favoured by studies of the proton structure function, F2. 

These predictions qualitatively describe the steep rise of the cross section with energy. 
At W = 250 GeV, the gluon density is being probed in these models [6] at x ~ 10~^, 
outside the range in which it is well constrained by global PDF analyses; the results are, 
therefore, sensitive to the PDF used. However, no discrimination between the gluon PDF 
can be made from the curves shown in Fig. 5 due to the large theoretical uncertainties 
from higher-twist contributions and the QCD scales due to missing higher-order terms. In 
addition, skewed parton distributions [51] are, as yet, relatively poorly constrained since 
the proton structure function F2 is not very sensitive to them. 

Also shown in Fig. 5 is the result (dot-dashed curve) of a calculation [52] based on a dipole 
model [53]. The J ftp wave-function was assumed to be Gaussian in both the transverse 
and longitudinal momenta of the quarks. The normalisation was fixed from the bo value 
reported in the next section. The W dependence of the model is in reasonable agreement 
with the present data. 

10 Differential cross-section da/dt and the Pomeron 
trajectory 

The differential cross-section da^yp^ji^p/ dt was calculated, in bins of W, separately for 
the muon and electron J /if) decay channels in the kinematic range -t < 1.8 GeV^ and 
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—t < 1.25 GeV^, respectively. The results are shown in Fig. 6 for four representative 
ranges of W. In each W bin, a fit of the form da/dt = da/dt\t=o ■ e''* was performed. 
For the muon sample, the fit was performed in the restricted range -t < 1.2GeV^ where 
the uncertainty resulting from the subtraction of proton-dissociative events is small. The 
results of the fits are given in Table 1. The muon and electron analyses give consistent 
results for da/dt\t=o and b, as shown in Figs. 6 and 7. 

The b slope increases with W and, in the geometrical picture of the interaction, is ap- 
proximately equal to that expected from the size of the proton [33], which suggests that 
the size of the J/ip is small compared to that of the proton. 

A value of a'jp was obtained by fitting the W dependence of b to the function b{W) = 
bo + Aa'jp ■ ln(iy/90 GeV), according to Eq. (2). The results for both the muon and elec- 
tron analyses, given in Table 2, are in good agreement. The systematic uncertainties were 
estimated by repeating the fit for each uncertainty not correlated in W and adding the de- 
viations from the nominal value in quadrature. The result of the combined measurement, 
shown as the line in Fig. 7, is 

bo = 4.15 ± 0.05{stat.)toil{syst.) GeV"^ 

a'jp = 0.116 ± 0.026(stat.)+°;°^°(syst.) GeV'^ 

The systematic uncertainties were computed from the combination of the muon and elec- 
tron analyses, taking into account the common systematic uncertainties. 

The Pomeron trajectory was determined directly by measuring the variation of the W 
dependence of the elastic cross section at fixed t, as parameterised in Eq. (1). This method 
is insensitive to the proton-dissociative background, since the latter was measured to be 
independent of W, as described in Section 6. In Fig. 8, the measurements of da/dt used 
in this determination of a]p{t) are presented as a function of W for fixed t; the line in 
each plot is the result of a fit of the form da/dt oc l^^ ['^iP(*)~^l. The resulting values 
of ajp{t), given in Table 3, are shown in Fig. 9, with the published HI results [15], as a 
function of t. They were fitted to the linear form ajpit) = a]p{0) + a'jpt. The separate 
fits from the muon and electron analyses are given in Table 2 and are in good agreement. 
The combined measurement gives 

aip(O) = 1.200 ± Om9{stat.)tofwi'Wst.) 

and 

a'jp = 0.115 ± OM8{stat.)tofilisyst.) GeV'^. 

The systematic uncertainties were computed from the combination of the muon and elec- 
tron analyses, taking into account the common systematic uncertainties. The result of 
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the fit was stable with respect to changes in the t range used for the fit. The slope, a'jp, 
measured with the present data, is not consistent with zero and therefore indicates a small 
shrinkage. The increase with W of the cross section and of the b slope, parameterised 
by aip{0) and a'jp, respectively, are in agreement with pQCD-based models [6,54]. The 
values of a]p{0) and 6 are compatible, after taking account of the measured value of a'jp. 
The soft-Pomeron trajectory ajp = 1.08 + 0.25 ■ t [11,55] is inconsistent with the present 
data. However, the contribution of the soft Pomeron plus a hard Pomeron as proposed 
by Donnachie and Landshoff [56] may well be able to describe the data. 



11 Conclusions 



The exclusive photoproduction of J/tp mesons has been studied at HERA with the ZEUS 
detector in the kinematic range 20 < W < 290 GeV using both the muon and the electron 
decay channels. The J/ip spin-density matrix elements, Tqq and rf^^ , have been measured; 
their values are consistent, within two standard deviations, with the hypothesis of s- 
channel helicity conservation. 

The 7p J/ijjp cross section exhibits a strong dependence on W, which can be parame- 
terised by a power- like dependence of the type W^, with S = 0.69±0.02{stat.)±0.03{syst.). 
This behaviour is described by pQCD-based models and can be understood as due to the 
increase of the gluon density in the proton for decreasing values of the parton fractional 
momentum. 

The differential cross-section da^p^j/^p/dt has been measured as a function of W for 
\t\ < 1.8 GeV^. It can be described by an exponential function in t, with a slope b = 
4.15 ± Om{stat.)tof8{syst.) GeV"^ at Vl^ = 90 GeV, which increases logarithmically with 
W. 

The parameters of the Pomeron trajectory, aip(O) and a'jp, have been determined from 
the W and t dependence of da^p^jj^p/ dt. The intercept is ajpiO) = 1.200 ± 0.009(stat.) 
toMtisyst.) and the slope is a'jp = 0.115 ± 0.018{stat.)tomlisyst-) GeV"^ These values 
are inconsistent with those expected from the exchange of a soft Pomeron. The data 
indicate that a'jp is different from zero but smaller by a factor of two than the value 
measured in soft hadronic interactions. 

Clearly therefore, the description of J/ ip production lies within the realm of perturbative 
QCD. A quantitative description comparable to the precision of the current data requires 
further theoretical progress. 
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Table 1: Measurements, in different ranges ofW, of the total J/tp photoproduction cross section, of the differential cross section 
extrapolated to t = and of the slope parameter b of the exponential t dependence. The first uncertainties are statistical 
and the second are systematic. Nobs is the number of events in the signal mass region, Nnon-res is the estimated non- 
resonant background, N^(2S) is the number of ipC^S) events in the J/ip mass region and A is the acceptance. The 
effective photon flux, $, is used to compute the 'jp cross section from the ep cross section. 
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J/ip ^ jjT jJ, 
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s 
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0.67 ±0.03 ±0.05 
30 < < 170 GeV 
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Eq. (2) 
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0.136 ±0.031+°:°°^ 
-t < 1.25 GeV2 


Fits to 
Eq. (1) 



Table 2: Measurements of 6, Bq, a'jp and aip(O) obtained separately from the 
muon and electron decay channels. The first uncertainties are statistical and the 
second are systematic. Where given, the third refers to the modelling of the proton- 
dissociative subtraction. The last column indicates how the values were determined. 
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Table 3: Values ofajp{t) obtained from fits to the function da/dt oc W^'^'^^^'^^ ^1. 
The first uncertainty is statistical and the second is the systematic. 
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Figure 1: Invariant-mass distributions of the ji^ jj, pairs in the different W 
regions: (a) 20 < < 30 GeV, (b) 30 < < 70 GeV, (a) 70 < W < 110 GeV, 
(d) 110 <W < 150 GeV and (e) 150 <W < 170 GeV . The histograms represent 
the LPAIR distributions of the non-resonant background. 
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Figure 2: Invariant-mass distributions of the e^e pairs in the different W 
regions: (a) 20 < W < 35 GeV, (b) 35 < W < 50 GeV, (a) 50 < < 90 GeV, 
(d) 90 < < 140 GeV, (e) 140 < < 200 GeV, (f) 200 < W < 260 GeV 
and (g) 260 <W < 290 GeV . The close-hatched histogram represents the LPAIR 
Monte Garlo distribution for the non-resonant background and the wide-hatched 
histogram that from G0MPT0N2. 
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Figure 3: T/ie acceptance- corrected decay angular distributions for exclusive J/ip 
photoproduction in the kinematic range 30 < W < 170 GeV and \t\ < 1 GeV'^. 
The non-resonant background has been subtracted. The results of both the fi'^fi~ 
and e"*"e~ decay channels are presented. The vertical bars indicate the statistical 
uncertainties only. The curves are the results of the fits to Eqs. (3) and (4), as 
described in the text. 
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Figure 4: The exclusive J ftp photoproduction cross section as a function of 
W for J/ip ^ fi~^fi~ and J/ip ^ e^e~ . The inner bars indicate the statistical 
uncertainties; the outer bars are the statistical and systematic uncertainties added 
in quadrature. Results from the HI [15], E4OI [44] O'nd E516 [45] experiments 
are also shown. The solid line is the result of a fit to the ZEUS data of the form 
a (X {W/90 GeVY and the dotted line is the extrapolation of the fit. 
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Figure 5: The exclusive J/ip photoproduction cross section as a function of W 
obtained in the two leptonic decay channels, J/tp fJ'~^f^~ (iiT-d J/ip ^ e^e~ . The 
inner bars indicate the statistical uncertainties, the outer bars are the statistical 
and systematic uncertainties added in quadrature. The experimental results are 
compared to the QCD predictions of MRT [5], using two different par ameterisations 
of the gluon PDF m the proton, MRST99 [48] (dotted curve) and CTEQ5M [47] 
(dashed curve). The solid curve shows the QCD prediction of FMS [6] using A = 4 
and the CTEQ4L [50] PDF. The dash-dotted curve displays the prediction [52] 
based on the colour-dipole model [53] with a double- Gaussian J/ip wave-function; 
this prediction was re-scaled to the h-slope measured in this paper. 
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Figure 6: The differential cross-section da^p^j/^p/dt for exclusive J/ip photopro- 
duction for representative bins of W and for the decay channels, J/ip fi'^n~ 
(squares) and J/ip — > e'^e" (points). The vertical bars indicate the statisti- 
cal uncertainties only. The full lines represent the results of a fit of the form 
da/dt = da / dt\t=Q ■ e^^ performed in the range —t < 1.2 GeV"^ for the muon channel 
and in the range —t < 1.25 GeV'^ for the electron channel. 
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Figure 7: Values of the slope, b, of the t distribution, plotted as a function 
of W . The line shows the result of a fit of the form b{W) = 6(90 GeV) + 4 ■ 
a'jpln{W/90 GeV). 
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Figure 8: The differential cross-section da^p-,j/^p/ dt as a function ofW at fixed 
t values. Only the statistical uncertainties are shown. The lines correspond to the 
results of fits of the form da^p^j/^p/ dt oc H/^['^Jp(*)~i] . 
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Figure 9: Pomeron trajectory as a function of t as obtained in the two lep- 
tonic decay channels, J/ip ^ A^"*"/^" Q^f^ J/i^ ~* The inner bars indicate 

the statistical uncertainties; the outer bars are the statistical and systematic uncer- 
tainties added in quadrature. The results from the HI Collaboration [15] are also 
shown. The solid and dotted lines are the result of linear fits to the ZEUS and HI 
data, respectively. The one standard deviation contour is indicated for the ZEUS 
(shaded area) and HI (dotted lines) measurements. The dashed line shows the DL 
soft-Pomeron trajectory [55]. 
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